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ABSTRACT 

An empirically motivated model is presented for accretion-dominated growth of the 
super massive black holes (SMBH) in galaxies, and the implications are studied for 
the evolution of the quasar population in the universe. We investigate the core as- 
pects of the quasar population, including space density evolution, evolution of the 
characteristic luminosity, plausible minimum masses of quasars, the mass function of 
SMBH and their formation epoch distribution. Our model suggests that the charac- 
teristic luminosity in the quasar luminosity function arises primarily as a consequence 
of a characteristic mass scale above which there is a systematic separation between the 
black hole and the halo merging rates. At lower mass scales, black hole merging closely 
tracks the merging of dark haloes. When combined with a declining efficiency of black 
hole formation with redshift, the model can reproduce the quasar luminosity function 
over a wide range of redshifts. The observed space density evolution of quasars is well 
described by formation rates of SMBH above ^ IQ^Mq. The inferred mass density 
of SMBH agrees with that found independently from estimates of the SMBH mass 
function derived empirically from the quasar luminosity function. 

Key words: cosmology: dark matter - galaxies: formation - galaxies: active - quasars: 
general 
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1 INTRODUCTION 

Super massive black holes consisting of at least tens of mil- 
lions of solar masses are now believed to be lurking in the 
nuclei of most galaxies. These dark objects are also the cen- 
tral engines of the most powerful sources of light in the 
universe: the quasars. Discerning the evolution history of 
SMBHs and quasars is a fundamental challenge for the the- 
ories of structure formation in the universe. Quasars and 
other forms of Active Galactic Nuclei (AGN) are thought 
to be powered by accretion on to the SMBHs. By virtue of 
their enormous brightness these objects are now observable 
up to sufficiently high redshifts, corresponding to the dawn 
of the era of galaxy and star formation in the universe. In 
the context of hierarchical structure formation, it is believed 
that mergers of galaxies play a crucial role in igniting the 
quasars. The merger induced gravitational instabilities fun- 
nel a huge amount of gas to the galactic cores. As a fraction 
of the kinetic energy of this gas is released in the form of 
electromagnetic radiation, luminosities close to the Edding- 
ton luminosity (~ lO"*' erg s'^ for atypical^ 10* M© SMBH 
masses) are produced. 
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Nuclear black holes probably trace the merger of their 
host galaxies and coalesce during mergers. If one envisages a 
simple scenario where a dark matter halo harbours a single 
galaxy and halo merging leads to galaxy merging along with 
their SMBHs, a coalescence dominated growth naturally 
leads to a linear relation between the black hole and the host 
halo masses (see e.g. Haehnelt, Natarajan & Rees 1998, here- 
after HNR98). The observed correlation between the black 
hole masses and t he bulge velocity disp ersion (Mbh — cr re- 
lation) of galaxies jPerrarese fc Merritti200Cl,l . however, sug- 
gests that the accreted component of mass far supersedes the 
coalesced one. The observed scaling is concomitant with a 
scenario where the energy feedback from the black hole plays 
a crucial role i n regulating the amount of matter accreted 
on to it llSilk fc Rees. .IQQSil . During the accretion process 
the mass of the SMBH builds up over some characteristic 
time-scale. The ratio of this characteristic time-scale to the 
Hubble time determines the duty cycle of a quasar. Accre- 
tion shuts down when the black hole is so massive that its 
energy feedback exceeds the binding energy of its fuel reser- 
voir. I ndeed, as pointed out by HNR98 and lWvithe fc Loebl 
i2003h . the observed scaling can be retrieved if the black hole 
were to supply the gas with an energy equal to its binding 
energy over the dynamical time of galaxies. These intricately 
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woven, yet neat set of ideas sketch an outline for simple mod- 
els of quasar and SMBH evolution. 

Several analytic and semi-analytic models based on 
mergers of dark haloes have been proposed to explain the 
high redshift evolution of quasar luminosity functions (see 
for e.g. Haehnelt & Rees 1993, hereafter HR93; HNR98; 
Haiman & Loeb 1998; Kauffmann & Haehnelt 2000, here- 
after KHOO; CavaUere & Vittorini 2000, hereafter CVOO; 
Hosokawa 2002; Wyithe & Loeb 2003, hereafter WL03; 
Hatziminaoglou et al. 2003). Simplified analytic models as- 
sume that all the cold baryons in a halo are locked in a single 
galaxy and galaxy mergers trace the halo mergers. More re- 
alistically, the progenitor galaxies in a newly formed halo 
sink towards the centre over a dynamical friction time-scale 
(Binney & Tremaine 1989; Lacey & Cole 1993) and merge to 
form a new baryonic core. However, as the dynamical fric- 
tion time-scale is generally much shorter than the Hubble 
time, the approximation that these mergers are instanta- 
neous, works rather well. Galactic mergers are assumed to 
ignite an AGN lasting for a characteristic time tqso- The re- 
sulting luminosities are estimated as Eddington luminosities 
corresponding to central black hole masses, which in turn 
are determined through their correlation with the circular 
velocities of the newly formed haloes. 

The uncertainty in the mass of cold gas available for 
accretion poses a key difficulty for the models. There are two 
central aspects of this uncertainty: the first is to ascertain 
how much col d gas is available in galaxies of a given mass at 
a given epoch jKauffmann fc Haehnehl200(J) and the second 
is to understand how effectively this gas is transported via 
the galaxy merg ers to become available for a ccretion by the 
forming SMBH JCavaliere fc Vittorinillioool) . 

The predictions of these models can be tested against 
the backdrop of a fairly well constrained observational data 
set. The most striking aspect of this data is that the quasar 
population peaks at z ~ 2.5 and exhibits a steep decline 
towards lower redshifts (Schmidt, Schneider & Gunn 1995; 
Warren, Hewett & Osmer 1994; Shaver et al. 1996; Boyle 
et al. 2000). At redshifts higher than ~ 2.5 the number 
count appears to decline again, however, the observations 
at these redshift could be marred by obscuration and ex- 
tinction effects. Simple analytical models mentioned above 
can well explain the high redshift observed luminosity func- 
tions {z > 2). They, however, systematically fail at lower 
redshifts {z < 2) and this epoch in the history of the uni- 
verse is somewhat opaque in our understanding of the evolu- 
tion of quasars. Even though the works of CVOO, HR93 and 
KHOO shed important light in this regard, no comprehen- 
sive scheme is yet available to systematically quantify the 
various aspects of quasar evolution at z < 2. 

The present work is primarily motivated by the low red- 
shift anomaly between observations and theoretical predic- 
tions. To bridge this gap we introduce a toy model cen- 
tred around two main ideas: the first pertains to the effect 
of declining black formation efficiency at low redshifts and 
the second involves a declining black hole formation rate in 
haloes above some characteristic mass M, . The first idea 
is der ived from a model investigated in iHaehnelt fc Ree'3 
il993h where it is argued that the conditions for SMBH 
formation are more favourable at higher redshifts as com- 
pared to low redshifts. The second idea is based on the 
existence of a break luminosity in the empirical quasar lu- 



minosity functions. We attribute this break luminosity to 
a systematic decoupling between halo and galaxy merger 
rates in haloes above a characteristic mass along the line 
suggested by CVOO that galaxy merging becomes rarer in 
sufficiently massive haloes, where individual galaxies re- 
tain their properties. Furthermo re, in a similar context, 
iHaiman. Ciotti fc Ostrikeil l|200^ have argued that a de- 
crease in the quasar luminosity density towards low redshifts 
arises due to a decline in the formation rates of spheroids. 
Here we propose a characteristic mass scale A/* for haloes 
which marks a drop in the formation rates of black holes 
and thus possibly spheroids as well. As the low redshift uni- 
verse is an era of group assemblage, our value of M, which 
corresponds to a group mass scale seems natural to explain 
the steep decline observed in the quasar luminosity density 
at low z. 

Our model entails a phenomenological treatment of the 
ideas described above. Notwithstanding the crudeness of the 
approach, it yields interesting results for the the quasar 
luminosity functions and their evolving number density, 
broadly consistent with the observations. It also provides 
a simple explanation for the evolution of the characteristic 
luminosity of quasars. Furthermore, the minimum SMBH 
masses for quasars suggested by the model (by comparison 
with the observed space density evolution of quasars) are 
> 10^ Mp), wh ich is consistent with the original estimates of 
ISohanI 1119831 . Finally the model can be linked quite simply 
with the evolution of cold gas masses in galaxies and star for- 
mation rates in the universe and thus yields estimates of the 
contributions to the ionising background at high redshifts, 
due to star formation and quasar activity. 

The structure of the paper is as follows: in section 2 we 
review some important aspects of analytical modelling of hi- 
erarchical structure formation which provides the essential 
tools for our model; section 3 discusses the relation between 
the black hole masses and their host haloes in a scenario 
with a declining efficiency of SMBH formation with redshift. 
Section 4 extends the argument of major mergers of haloes 
to compute the SMBH formation rates and the luminosity 
functions of the ensuing quasars. In section 5 we compare 
the SMBH formation rates from our model against those in- 
ferred from the empirical double power-law fit to the quasar 
luminosity functions. In section 6 we estimate the mass func- 
tions of SMBH and the comoving cosmological mass density 
in the SMBH at various redshifts. Section 7 describes how 
to obtain simple estimates of the cold gas density and star 
formation rate density in the universe; section 8 extrapolates 
these star formation rates to high redshift to compute the 
ionising background. These are then compared with a simi- 
lar extrapolation to high redshifts, of the quasar luminosity 
density. We summarise and discuss our results in section 9. 



2 ELEMENTS FROM EXTENDED 
PRESS-SCHECHTER THEORY 

Our approach largely takes advantage of the extended Press- 
Schechter (EPS) formalism which provides essential tools for 
computing the merger probabilities and rates of bound ob- 
jects (Bond et al. 1991; Bower 1991; Lacey fc Cole 1993, 
hereafter LC93). In this scenario the initial density fluctua- 
tions in the universe are assumed to be Gaussian in nature. 
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The growth of these fluctuations proceeds hnearly initially 
and is analysed through the linearised ideal fluid equations 
(e.g. Silk & Bouwens 2001). The fluctuations subsequently 
turn non-linear and eventually collapse to form bound ob- 
jects as the density reaches a critical threshold. A compar- 
ison of EP S results with the re sults of N-body simulations 
as given in lLacev fc Cold il994) . demonstrates that the pre- 
scription works rather well. 

The linear growth of fluctuations is described through 
the linear growth factor D{t). For a critical cosmology this 
factor is simply equal to the scale factor a for that epoch. For 
the ACDM cosmology used in this paper {Q,^ = 0.3, = 
0.7, = 0.02//t^ as = 0.9,/t = 0.7) . we compute it as 
given in lNavarro. Frenk fc White! 1^9^. The mass M of an 
inhomogeneity is represented through a sharp k-space fll- 
ter which is used to smooth the matter density held at this 
mass scale. As suggested in LC93, M can be related to ks 
through a qualitative mapping M = Gn^pokJ'^- Here po is 
the mean background density of matter in the universe. Us- 
in g the dark matter power-spectrum P{k) parametrized as 
in lEfstathiou. Bond fc White! (!l992h . the variance for mass 
scale M can then be simply calculated as 



a (M) = / P{k) 4tt k^ dk. 
Jo 



(1) 



The mass function of haloes is now written in terms of 
(j^(M) as 



NpsiM,t) 



2 po 



TT M D(t) 



1 d a 



exp 



2D2(t)cr2 



(2) 

where uj is the critical density for collapse. 

The formation rates, major merger rates and survival 
probabilities have been discussed by several authors (e.g. 
Sasaki 1994; Kitayama & Suto 1996, hereafter KS96). In this 
context, the excursion set approach of Bond ct al. (1991i) 
(also LC93) yields two important conditional probabilities. 
The first is the conditional probability that a halo of given 
mass M2 at time t2 , had a progenitor of mass Mi at earlier 
time ti. This is given as 



Pl{Sl,UJi\S2,lJJ2)dSl = 



1 



- UJ2 



2^ (5l - S2)3/2 

{lui — UJ2Y 



xexp 



2(Si -S2 



dSi 



(3) 



where U}\ is the critical density at time t\ and C1J2 at time t2- 
Also, 5i is the variance corresponding to mass M\ and S2 is 
the variance corresponding to M2. The relation between the 
critical density w(t) at any time t is given as tj(t) = ujq/ D{t), 
where loq is the present day critical density (-^ 1.686 for a 
SCDM cosmology). 

The second is the probability that a halo of mass Mi at 
time t\ , will result in a larger halo of mass M2 at later time 
t2 can be written as 



P2{S2,i^2\Sl,UJl)AS2 



1 LJ2(ti'l — iJJ2) 



2-K 



xexp 



LJl 
(S2a;i 



Si 



1 3/2 



S2{Sl-S2) 
SlWz)'" 



2SiS2(Si — S2 



dS2 



(4) 



Through appropriate transformations the conditional prob- 
abilities can also be written in terms of halo masses. 



2.1 Merger rates, survival probabilities and 
formation epoch distributions 

The EPS formalism lacks a self sufficient way to compute 
the formation rates of haloes. This is because effectively all 
haloes at any epoch are continuously accreting mass and in 
this sense are newly created. There is no clear distinction be- 
tween the objects accreting less mass and the objects which 
accrete substantial mass. For computing the major merger 
rates, therefore, a criterion has to be introduced by hand. 
One has to assume that if none of the progenitors of a newly 
formed object had mass greater than half (or some other 
fraction) of the object mass, the event is a major merger. 
Alternatively Sasaki (1994) proposed to derive the forma- 
tion rates by assuming that the destruction probability of 
objects (i.e. the probability that they get incorporated into 
higher mass objects) has no characteristic mass scale. This 
rate is then simply the positive term in the time derivative 
of the PS mass function. In Sasaki's prescription the survival 
probability of haloes from one redshift to another is simply 
the ratio of the growth factors at these redshift. The prod- 
uct of formation rate and the survival probability defines 
the formation epoch distribution. Interestingly for a criti- 
cal cosmology, Sasaki's formation epoch distribution can be 
integrated analytically over all redshifts to yield the exact 
mass function of haloes. 

In the present work we compute the major merger rates 
based on the prescriptions described in LC93 and KS96.The 
probability that a given halo of mass M2 at time t2 had, 
at some earlier time ti, a progenitor with mass below some 
threshold mass Mc(< M2) is computed by integrating equa- 
tion |3] over Si from S{Mc) to 00. If the difference between 
the times ti and t2 is substantial, then Mi can reach the 
value M2 through either a single merger occurring anytime 
between ti and t2 or it can do so through more than one 
merger /accretion episodes between the two times. However, 
if the difference between the two times is very small then the 
entire change in mass (AM = M2 — Mi) can be treated as 
arising from a single merging event. Letting Aui — ui — UJ2 
we thus write 



P{Mi < M,,ti\M2,t2) = 



= erf 



Pi dM, 



(5) 



^2(5. - 52) 



where Sc is the variance corresponding to mass Mc- For 
small Au! taking the series expansion in equation |S] up to 
first order only, we compute the major merger rates of dark 
haloes as 



RMM{M,t) 



Utt 



Nps{M,t), 



TT ^S{M/2) - S{M) 

(6) 

where we have chosen Mc = M/2. Note that when Sc tends 
to 5", the right hand side of equation |^ tends to one, if Au; 
is small but finite. This implies that at any given epoch, all 
the existing objects continuously accrete mass. 

As pointed out in LC93, when the halo mass for a tra- 
jectory falls to a smaller value at an earlier time, it only 
means that one of the parent haloes had this smaller mass. 
As a result, it is always possible that the largest coeval pro- 
genitor already has more than half the mass of the final halo 
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Figure 1. The formation epocli distribution of DM haloes of 
different masses as given by equation 151 



at this epoch. Such an event cannot be characterised as a 
major merger as we define it. The approach of KS does not 
account for this possibihty. Ahhough this marginally alters 
the rates for the range of halo masses we are concerned with, 
it is nevertheless interesting to see how one can systemati- 
cally account for this possibility and therefore, we propose 
a solution to this problem in appendix A. 

Equation |1| can be used to calculate the survival prob- 
abilities as discussed in LC93 and KS96. These are based 
on the assumption that an object retains its identity up to 
the time it doubles its mass. The probability that a halo 
of mass M at present time ti, will not undergo such a 
merger /accretion, before time t2, which would result in its 
final mass being greater than 2M, is given as 

Psurv{MMM)= P2{M2,t2\M,ti)dM2. (7) 

J M 

As discussed previously, the survival probability of 
haloes multiplied with their merger rates, yields their for- 
mation epoch distribution. This distribution tells us that in 
a given population of mass AI objects at time t, what frac- 
tion were formed at an earlier time tf and have retained their 
identity since then. Thus the formation epoch distribution 
in terms of formation time tf can be written as 



F{M,tf t)dMdtf 



RMM{M,tf) X Psurv{M,tf t) 

X AM dtf. (8) 

We show in FigQthe formation epoch distribution of haloes 
of different masses as a function of redshift, and in Fig|5|the 
formation rates of objects with masses above lO", 10^^ and 
10^^ Mq respectively. It can be seen from Fig|5|that merger 
rates of haloes alone do not produce a significant evolution 
in the number density at low z, and therefore a naive model 
directly linking quasar triggering to merger rates of haloes 
would fail to match the steep decline which is observed in 



Figure 2. The evolution of cumulative formation rates of dark 
matter haloes with masses above Mmin- These rates (particularly 
for Mmin = 10^^ Mq, the typical quasar host halo mass) do not 
produce as steep a decline at low redshifts as that measured for 
quasar number counts. 



their number density. Hence, we subsequently discuss the ad- 
ditional inputs that need to be incorporated into the model 
in order to reproduce the observed evolution of the quasar 
population. 



2.2 Halo circular velocities 

For a nearly isothermal halo, the mapping between the halo 
circular velocity and its mass can bo obtained using the re- 
lations given in Somorvillo & Primack ( 1999) 



GM f^A 2 2 



(9) 



where z is the collapse redshift and r^ir is the virial 
radius, be computed by noting that Ad = 

47r/3 r^i^(z) pvir{z) where the virial density pvir{z) is given 
as 



pvir{z) = Ac(2:)p°(l + Z) 



3 ^0 

Q.{z) 



(10) 



Here, Qq and Q,{z) are the matter densities of the present 
day and redshift z universe respectively; ~ 2.8 x 
10^^ /i^ MqMpc"^ is the present critical density of the uni- 
verse and a fitting expression for Ac{z) i s given as Ac(z) = 
IStt^ -I- 82a: — 38x^, where x = Qmiz) — 1 jSrvan fc NormanI 
T99^. Thus the expression for virial radius in terms of for- 
mation redshift z and halo mass Af is given as 



2.14 X lO'^AI 



1/3 



A,{z) 



nm{z) 187r2 



-1/3 



il + z)~^ Mpc. 

(11) 
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Figure 3. Phenomenological redshift evolution of the formation 
efficiency of black holes (see text for the exact function). The 
characteristic redshift of decline is Zt = I. A steep decline is 
obtained around 2 = where values become as low as ten percent 
of the maximum value. 



These relations together enable us to write the halo circular 
velocity as 



1.38 X 10"^M 



2 ,,1/3 



1/6 



(12) 



3 SUPER MASSIVE BLACK HOLE MASSES 

In a self-regulated growth of SMBH masses, the correla- 
tion between the galaxy circular velocities (uc) and SMBH 
masses {Mbh) is written as (.Silk fc Rees..l99&'l 



Mbh oc «c 



(13) 



where the constant of proportionality depends on the frac- 
tion of cold gas fcoid available for accretion by the coalesced 
seed black hole and weakly on the presence of an accretion 
disc (HNR98; Burkert & Silk 2001). Such a scaling is corrob- 
orat ed by the investigatio ns of p^'errarcsc & Mcrritt (2000) 
and lGebhardt et alj J2000h . 

Incorporating fcou in an analytic model is not simple. 
Firstly, we do not know how much cold gas exists in haloes 
of different masses at different epochs and secondly, even 
if we knew the mass of cold gas, it is impossible to ana- 
lytically compute the contributions coming from different 
mass progenitors during a merger. For the sake of simplic- 
ity, therefore, we get around this problem by adopting an 
empirically motivated decline in black hole formation effi- 
ciency with redshift. This can be attributed to a general 
decline in the amount of cold gas in the galaxies at low red- 
shifts. Such an approach has been co nsidered previously, for 
example, bv lHaehnelt fc ReesI (Il993l) . 



We model the declining formation efficiency of SMBH 
masses through a function e{z) such that Mbh ~ ^{z) 
and consider a phenomenological form of e[z) as 



e{z) = eo X 



c — exp 



(14) 



Here eo and c are constants and is a characteristic redshift. 
We choose a value a, ~ 1 based on the observed peak in 
the star formation rate (SFR) density in the universe at 
redshifts between 1 and 2 (see e.g. Ascasibar et al. 2002; Silk 
& Rees 1998) . This is because a peak in the SFR density is a 
general indicator of a gas rich environment in galaxies. It is, 
however, still a matter of debate whether SFR density peaks 
at z ~ 1 — 2 or reaches a broad maximum there. To obtain 
a steep decline we choose c = 1.1 such that the minimum 
value of e(z) is roughly 10 percent of the maximum value 
1.1 X eo ■ This latter is chosen to match the observed Mbh^o 
correlation (Ferrarcsc 2002) 



Mbh = 1.66 x 10** 



Me 



(15) 



200 km s-i 

where we adopt a = 5 consistent with a feedback regu- 
lated growth of SMBH masses. This is slightly different from 
the value a — 4.58 ± 0.52 inferred from SMBH mass mea- 
surements, the difference possibly having its origin in the 
observed non-linearity of Vc — a relation (Ferrarese 2002; 
WL03). For simplicity, we consider the linearly fitted rela- 
tion for ellipticals a « 0.65 Vc to obtain eo = 6 x 10""" (such 
that Mbh ~ eo at high redshifts). 

Figure 13 shows a plot of efficiency versus the formation 
redshift, i.e. how eauation ll4l behaves. We can now write the 
black hole masses in terms of formation redshift z and halo 
masses M as 



Me 



e{z) X 5 X 10-^° — 



M 

Mq 



5/3 



Ae(z) 



3/6 



(1 + 2 



,5/2 



MfT 



(16) 



4 FORMATION RATES OF SMBH AND 
QUASAR LUMINOSITY FUNCTIONS 

We now use the Mbh — Mhaio mapping (equation 1161 to 
transform the halo major merger rates into the black hole 
formation rates in terms of the formation redshift z as 



B)bh{MbH,z) = RMM[Mhalo{MBH,z),z\ 



dM 

dMBH 



{Mbh,z). 
(17) 

The superscript 1 in B^bu is to differentiate it from a mod- 
ified rate Rbh, which takes into account a decline in black 
hole formation rates in haloes above a characteristic mass 
Mf This is achieved through a phenomenological function 

Ff_,{Mhaio) such that 

Rbh {Mbh, z) = Rbh {Mbh, z) x F4Mhaio{MBH , z)]. 

(18) 

For the form of Ff^{Mhaio), we found that a simple expo- 
nential function exp(— M/M»), with a cutoff at M« is too 
steep to be consistent with the observed luminosity func- 
tions. We therefore adopt a function with a following be- 
haviour: F{M) 1 for M « M, and F{M) {M^/MY 
for M » M, , where is a positive number. Instead of a 
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Figure 4. This figure shows the second of our two phenomenolog- 
ical functions. We use this function to introduce a drop in the for- 
mation rates of black hole in haloes above the characteristic mass 
scale Mt , the value of which has been set by comparison with the 
quasar luminosity function l'figure l2Tl to about M* = 
The redshift evolution of the luminosity corresponding to this 
characteristic mass is plotted against the empirical characteristic 
luminosity of the double power law fit, in figlTI 



sharp exponential cut-oflt, such a function has softer 'power- 
law' drop at the high mass end: 



(19) 



1 + {M/M,y' 

Ffj,{AI) for values of /i = 1 and ^ = 2 is plotted in figure 3] 
along with the exponential function F{M) = exp(— M/M*). 

It is important to mention at this point that in mas- 
sive haloes (M > M*), rarer galaxy mergers could possibly 
alter the Mbh — Mhaio mapping. Furthermore, the AGN 
could enter the sub-Eddington regime due to lack of gas 
mass transported through rare galaxy mergers. However, to 
keep our model simple we assume that even for ineffective 
galaxy mergers, the MBH — Mhaio scaling is as given in equa- 
tion^] If a merger-ignited black hole shines at the Edding- 
ton luminosity Lf,ai/LEdd = 1, we can write the i?— band 
luminosities of AGNs as 



Lb ^ Jb X 1.3 X 10^ 



Me 



J s" 



(20) 



where /s is the fraction of bolometric luminosity emitted 
in the B-band. Realistically, as the accretion proceeds, the 
black hole mass builds up over the characteristic accretion 
time. This process is accompanied by a parallel increase in 
the quasar luminosity which determines its light curve. How- 
ever, a simplification is possible if one considers that the 
black hole mass instantaneously increases to post merger 
values given by equation 1161 The black hole then shines at 
the Eddington luminosity, corresponding to this new mass, 
over characteristic time-scale tqso which lies in the range 
tqso/tHubbie ~ 10~^ to 10~^. Thus usiug the Lb — Mbh 



Figure 5. The model lifetime of quasars given as tqso = k tdyn, 
where t^iyn is the halo dynamical time and the constant of pro- 
portionality is about k 0.035 (e.g. WL03) 



mapping and the fact that tqso <C tnubbie we can simply 
write the quasar luminosity function as 

^{LB,z)~tqsoy.RBH[MBH{LB),z\y.'^^^^^{LB,z). (21) 

The quasar lifetime, tqso is not a very well constrained 
quantity. Two relevant time-scales for quasars are the 
Salpeter (or e-folding) time-scale and the dynamical time- 
scale of galaxies. The Salpeter time-scale is given by tsai = 
4 X IQ^ {trad/Q-^) where erad is the radiative efficiency of 
black holes defined through the relation L^,ol — eradMc? , c 
is the speed of fight and M is the accretion rate of mat- 
ter on to the black hole. On the other hand it has been 
noted, for example, in HNR98 and WL03, that the observed 
Mbh — Vc relation naturally arises if the black hole injects 
the gas with an energy equal to its binding energy, over a 
local dynamical time-scale of the galaxy. In this respect, the 
local dynamical time of a galaxy is also a natural time-scale. 
In the present model, following WL03, we adopt a quasar 
lifetime proportional to the local dynamical time of galaxies 
in newly formed haloes. The local dynamical time of galaxies 
being themselves proportional to the halo dynamical time, 
we get tqso oc tdyn, where tdyn is the halo dynamical time 
and so we write 

-1/2 

(l + 2)-^/Vr, (22) 



tqso — kxl.5x 10' 



A,(2 



0.035 (e.g. 



where the constant of proportionality is k 
WL03). 

In figure |S| we plot the luminosity functions as given 
by equation 1211 (solid curves). The parameters are ^ — 2, 
M, = io"-25 Mq and /s = 0.06. The short- dashed curve is 
the double power law fit as given in^ej 1^9^. We transform 
the data and the double power-law fit to a A cosmology by 
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Log(LB/Lo) 



Figure 6. Model luminosity functions (solid curve) plotted against the observed luminosity functions (filled circle) of quasars. The data 
for 2 < 4 is fromlPei 1 1995) and for z Ri 4.4 is from Kennefick, Diorgovski & Carvalhc 1 1995). The model p arameters are eo = 6 X 10~^, 
11 = 2, Af. = 1013-25;v/q and 

iqso ^ 0.035 tfiyj^. The short-dashed curves show the double power law fit of iPeJ 119951) . The long-dashed 
vertical line is the model break luminosity corresponding to M» and the short-dashed vertical line is the characteristic luminosity in the 
double power-law fit. 



noting that oc dV~^dJ^^ and L oc di, where dV is 

the volume element and is the luminosity distance (e.g. 
Hosokawa 2002). The double power-law is given as 



(23) 



where log($./Gpc-=') = 2.95 ± 0.15, Pi = 1.64 ± 0.18 and 
Ph ~ 3.52 ± 0.11. The characteristic luminosity Lz has a 



Gaussian form 

=L.(l-f z)-'i+"^exp 



2a! 



(24) 



where log(i*/i0) = 13.03±0.10, = 2.75±0.05, a = -0.5 
and CT« = 0.93 ± 0.03. In figure El the vertical short-dashed 
line is the characteristic luminosity given by eauation l24l and 
the vertical long-dashed line shows the break luminosity for 
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Figure 7. The evolution of break luminosity in our model (solid 
line) compared with the characteristic luminosity in the double 
power law luminosity function fit (dashed line). Also plotted is 
the break luminosity in a scenario where e = eo is fixed (dotted 



line). The characteristic mass is M« 
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Figure 8. This figure shows the redshift evolution of formation 
rates of black holes with masses 10^,10^,10** and 10^ Mq respec- 
tively. The dashed lines are from equation 1251 and the solid line 
are from equation 1181 



Figure 9. Cumulative formation rates of black holes with masses 
above 10^, 10^ and 10^ Mq respectively. 



The characteristic luminosity and the break luminosity fol- 
low each other closely up to z ~ 2.5. This is especially clear 
in figure [7| where we additionally show (dotted curve) the 
break luminosity corresponding to A/* in a fixed black hole 
formation efficiency scenario. 

It has been suggested in WL03 that rather than be- 
ing an intrinsic feature of the quasars themselves, the break 
in the luminosity function is present only at low redshifts. 
Our model suggests otherwise. We find that the break in the 
luminosity function exists as a consequence of the charac- 
teristic mass and is present over a wide range of redshifts. 
Furthermore, the double power law in the form given above 
under-predicts the quasar luminosity function at 2 = 4.4. 
On the other hand, inserting the break luminosity given by 
our model into the double power law description of the data 
over-predicts the quasar luminosity function. A possible im- 
plication of this is that perhaps the black hole formation 
efficiency declines beyond 2 ~ 4. 

The tuning of parameters for our model luminosity 
function is admittedly crude. Also, it can be seen from equa- 
tions^] and I^U] that there exists a degeneracy between the 
parameters eo and fs ■ Higher values of eo require systemat- 
ically lower values of js- Nevertheless our fiducial numbers 
are in broad agreement with the values obtained by other 
investigations, and we emphasize once again that the model 
presented here only attempts to discern the broad trends of 
SMBH and quasar population evolution. 



our model, corresponding to Af» = lO'^^ '^^Af© haloes. It can 
be seen from the figure that the model can systematically 
account for the observed data over a wide range of redshifts. 
At 2 < 1, where the break luminosity of the model is slightly 
greater than the empirically fitted characteristic luminos- 
ity, we find that the bright end is slightly over predicted. 



5 FORMATION RATES OF SMBH FROM THE 
OBSERVED LUMINOSITY FUNCTIONS 

Due to their emission of gravitational waves, the merger 
rates of SMBHs will be detected directly by instruments such 
as the Laser Interferometer Space Antenna (LISA). In our 
model, SMBH merger/formation rates are given by equation 
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Figure 10. The space density evolution of quasars with minimum black hole mass Mmin in the model, plotted against the observational 
data. Filled circles describe the observed space density evolution of quasars with absolute B-band magnitudes Mb < —24 and open 
triangles describe the observed space density evolution of quasars with absolute B-band magnitudes Mg < —26. The data points are 
plotte d as in KaufFmann & Haehnclt 2000 (sec references therein). Also shown are the best fit to the high redshift data from lFan et alJ 
l200l|) (dashed lines) and .Schmidt. Schneider fc Gunn. L1995.) (dotted lines). The model parameters are same as in figure 151 



1181 One can also use the double power-law fit to the lumi- 
nosity function data to estimate these rates as a function of 
SMBH mass. These can be written as 

R°i!l{MBH,z) = ^ X <!>obs[LB{MBH,z),z] 

In figure IHl we plot SMBH formation rates for different 
masses using these two previous methods. The solid curves 
are computed from eauation ll8l i.e. our fiducial model, and 
the short-dashed curves are results from equation 1251 The 
different masses considered are labelled in the figure. It can 
be seen that there is a significant discrepancy at high red- 
shifts between the double-power law predictions and our 
model's. If, as discussed previously (section 4) , the efficiency 
e declines beyond z ~ 4, a steeper decline in the SMBH for- 
mation rates at high redshifts would be obtained. This would 
of course bring the two estimates in better agreement. In 
other words, this would mean that our model only gives the 
maximum limiting value of SMBH formation rates at these 
redshifts. 

In any case, it should be noted from the figure that 
the formation rates corresponding to 10^ and 10^ Mq black 
holes are predicted to peak at very high redshifts {z > 10) 
and have greater values (at least a couple of orders of mag- 
nitude larger) than their more massive counterparts. There- 
fore, their detection (or lack of detection) by LISA would 
yield strong constraints on our understanding of the forma- 
tion history of SMBH and their plausible minimum masses. 

Figure|5|shows the cumulative formation rates of SMBH 



computed by integrating eguation llSl from a minimum black 
hole mass Mmin to a tentative upper limit of ^ IO'^^Mq. 
Note that the decline in the cumulative formation rates of 
SMBH towards low redshifts is much steeper as compared to 
the drop in the cumulative formation rates of haloes (figure 

Hi. 

The number density evolution of the quasars with 
SMBH masses above some minimum mass Mmin can now 
be written as 

NgsoiMBH > Mmi„,z) = / dMBH X {tqsoiz) 

X Rbh(Mbh,z) }. (26) 

Figure irni shows the the number density evolution of quasars 
with minimum black hole masses as labelled in the fig- 
ure. Also plotted is the observed data for the space den- 
sity evolution of quasars with minimum B-band magnitudes 
Mb < -24 (filled circles) and Mb < -26 (open triangles) 
up to a redshift of 10. As can be seen, the figure sketches 
a consistent picture with the original est imate o f minim um 
quasar black hole masses of ~ 10* Mq bv lSoltanI (Il982h . 



6 THE MASS FUNCTION AND 

COSMOLOGICAL MASS DENSITY OF 
SMBHS 

The formation rates of SMBHs can be used to compute their 
formation epoch distribution (see section 2.1). These can 
then be integrated over redshift to yield the SMBH mass 
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Figure 11. The mass funct ion of SMBH s. Ope n squares show 
the local estimates given in ISalucci et al.l il999l) . The limiting 
2 = 2.5 mass functions from our model have been plotted as solid 
line (using Ng^) and long-dashed line (using N^jj). We also 
show the upper limit mass function at 2 = 1 (using Ng^) and at 
z = 6, where both limiting mass functions (Ng^ and Ngg) are 
very similar. 



T 




z (Redshift) 

Figure 12. Comoving cosmological mass density of SMBHs by 
using the SMBH mass function Ngjj (dotted lines) and Ngg 
(solid lines). Minimum SMBH masses NlYnin ctre (top to bottom) 
10*, 3 X 10* and 10^ Mq. The two curves (solid and dotted) can 
be considered as (upper and lower) limits on densities for each 



function. The formation epoch distribution of SMBHs can 
be written as 

FiH{MBH,zf z) = Rbh[Mbh,z;) X -^^(zj) (27) 

y<Psurv[Mhalo{MBH,Zf),Zf —* z] . 

However, there are two problems with this equation. The 
first is that the survival probabilities of SMBHs are only 
crudely correlated to the halo survival probabilities: the de- 
clining merger rates of SMBH in haloes with mass above 
Mt suggest that the survival probability may be greater for 
SMBHs in these haloes. The second problem is more fun- 
damental and lies in the fact that equation 1161 is not only 
a function of halo mass but also of formation redshift. This 
can be understood in the following way: masses of SMBH for 
a given mass halo will be higher at high redshifts (this holds 
even if the formation efficiency is taken as constant such 
that e{z) = eo), so that when we compute the destruction 
probability Pdest of objects (note that Psurv = 1 — Pdest), we 
authorise un-realistic situations where a halo with a higher 
mass black hole goes on to form a lower mass black hole. 
Hence we underestimate the survival probability. Equation 
1271 therefore, only gives a lower limit on the number counts 
of objects of different masses. An upper limit can be set by 
taking Paurv = 1. This amounts to computing mass func- 
tions by simply integrating the formation rates up to the 
redshift of observation. Such an approach has also been con- 
sidered in Hachnclt & Rcos (1993). Thus we write 

FEh{Mbh, zf z) = Rbh{Mbh, Zf) X -^(zf), zj ^ z. 

azf 

(28) 

We denote the corresponding SMBH mass functions as 
Nbh{Mbh,z) and Nbh{Mbh,z) respectively. These are 
computed as 

Nbh{Mbh, z) = r dzf F^if{MBH,Zf ^ z) . (29) 

•J OCi 

In figure 1111 we plot the SMBH mass functions for sev- 
eral redshifts and compare them with the estimates of local 
SMBH mass function as given in ISalucci et all |l99i). At 
redshifts greater than z ~ 3 there appears to be no signifi- 
cant difference in the two limiting mass functions {Nbh a^nd 
Nbh)- We can now write the comoving cosmological SMBH 
mass density as 

pi'^iM > Mrmn,z)= / dM B H {M B H X N^'j^ {M B H , z)} 

(30) 

Figure m\ shows the redshift evolution of cosmological 
mass density of the SMBH with minimum masses 10*, 3x10* 
and 10® Mq (top to bottom) respectively. The solid lines 
are computed using Ngg and dotted lines are computed 
using Nbh- It can be seen that mass function Nbh gives 
an unrealistic decline in the mass density of SMBH which 
should increase monotonically. Once again, the reason is that 
the survival probability of the SMBHs are underestimated 
and thus the number counts are underestimated as well'^. 



^ We point out tha t in a recent paper 
iBromlev. Somerville &: FabianI l2003l) have found a similar 
non-monotonic evolution of SMBH mass density, computed by 
using the M — Vc relation combined with the Press-Schechter 
mass function. 
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Figure 13. The contribution to the present day 2 = relic black holes from merger activity at higher redshifts. Different line-styles show 
results of the model for different black hole masses whose values are indicated in the figure. Each pair of curve in the same line-style 
indicates estimates using Fgjj (lower curve) and Fgjj (upper curve) as detailed in the text. It can be seen that at 2 = 0, a significant 
contribution comes from 2 > 1. 



The density computed using the mass function N§u thus 
sets an upper limit to the density of SMBHs for the three 
different minimum masses. 



In section 3 we have presented a scenario with a mono- 
tonically declining efficiency of black hole formation with 
decreasing redshift, with Mbh — Vc relatio n written as 
Mbh = t'c- However, the recent results of lShields et alJ 
i2003h point towards a redshift independence of Mbh — cr 
relation (and therefore, by extension, Mbh — Vc relation). 
They find no discernible change in the relation up to redshift 
of 2 « 3. Thus they suggest that black holes are well formed 
by « « 3. In order to reconcile the redshift dependence of 
our model Mbh — Vc scaling with observed redshift inde- 
pendence of this relation, we present an argument based on 
the formation epoch distribution of the SMBH. In figure 
we plot the formation epoch distribution of SMBH of dif- 
ferent masses at different redshifts of observation Zobs- It is 
evident from the figure that most relic SMBH between red- 
shifts z — and z — 1 were formed at high redshifts 2 > 1. 
As per our model, at redshifts z > 1, the function e{z) does 
not vary significantly as z « 1 is a characteristic redshift 
of decline for e(z). Appreciable decline in e{z) occurs only 
at a < 1. Hence, our adopted redshift varying Mbh — Vc 
relation does not pose a problem with the observed redshift 
independence of the relation. Furthermore, it can be seen 
from figure lT^ that using Fbh our model predicts that most 
relic SMBH between z = and z = 1, with masses above 
10* Mq, are formed between z ~ 2 and 3. Whereas, using 
Fbh we predict that most SMBH were formed at z ~ 1. 



7 ESTIMATES OF STAR FORMATION RATE 
AND NEUTRAL GAS DENSITY IN THE 
UNIVERSE 

In our model, the declining black hole formation efficiency 
is attributed to a corresponding decline in the amount of 
cold gas mass available in galaxies. To estimate this cold 
gas mass, we take the fraction of cold gas fgas[z) (which is 
different from the post merger fraction of cold gas available 
for accretion - fcoid) in galaxies as proportional to e{z), i.e. 
fgas{z) oc e(z). For simplicity we assume that this fraction 
is independent of halo masses. We fix the constant of pro- 
portionality to match the amount of cold gas observed at 
high redshift in Damped-Ly-Q systems . The cooling of gas 
in massive haloes may be ineffective (iB enson et alJli'oOSh 
and therefore massive haloes may be deficient in cold gas. 
We, therefore, tentatively adopt an upper limit to the halo 
masses which are rich in cold gas as cor- 
responding to cluster scales. 

Using these simple prescriptions, we can write the mass 
of cold gas as Mgas{M, z) oc e(z)M, where M is the halo 
mass. The comoving cosmological mass density of cold gas 
can then be written as 

Pgas{z)= dM {Mgas{M,z)xNps{M,z)}. (31) 

We take Mmin as the minimum halo mass cor- 
responding to the minimum halo circular velocity Vm 
for effect ively reta iinng gas in th e halo potential well 
( Somcrvi lle. Primack fc Fabe 3|2003). To simplify our com- 
putation we calculate Mmin corresponding to Vmin using 
equation 1121 bv setting the redshift of formation as the red- 
shift of observation itself. This approach has been adopted 
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Figure 14. The evolution of comoving cosmological mass den- 
sity of cold gas in galaxies as discussed in section 7. The data 
points are as given in lPeroux et alj (2001), and the different line- 
styles correspond to estimates truncated at different low circular 
velocities indicated on the figure (see text for details). 



previously, for example, in iRobinson fc Silkl i200Cl) . where 
they have discussed the error introduced by this simplifi- 
cation. Alternatively one can use the formation epoch dis- 
tribution and convert it from masses to circular velocities 
and then integrate over all redshifts up to the redshifts of 
observation to get the circular velocity functions. 

Dividing pgas by the critical density pc of matter in the 
universe we obtain flgas- This is plotted in figure [Tl] against 
the o bservations of Damped-Ly-a absorbers (jPero^^^^i] 
I2OOII) . As previously mentioned, we have normalised ilgaa 
to the highest redshift data point in figure fTH which yields 
a constant of proportionahty of the order of 10~^ for the 
fgas{z) — e{z) relation, with the exact value depending on 
the lower circular velocity cut-off used in the computation 

of ^gas • 

Further, the star formation rates can be estimated by 
adopting an SFR recipe as (e.g. Sommerville & Primack 
1999) 



M, 



Mgas 



(32) 



where e, is the efficiency of star formation and Mga3 is 
the cold gas mass available for star formation as described 
above. We can thus write the SFR density in the universe 
as 



Psfr{z) = 



dM M. X Nps{M, 



(33) 



Retaining the same normalisation as for the cold gas 
density, we choose tt = 10* years and an efficiency of 5 
percent i.e. e* = 0.05. Mmin and Mmax are as discussed for 



Figure 15. The evolution if SFR density in the universe. The 
different curves are results of the model with different minimum 
circular velocities using the same line-style code as in figure [T^ 
The long-dashed curve is from th e emp irical SFR density evo- 
lution given in lMadau &: P ozzetti 12000). The data points with 
error bars are as compiled iu iAscasibar et al.. 1.2002.) . 



the case of cold gas density. Figure 1151 shows the evolution 
of the SFR density in the universe. 

It can be seen from the figure that our crude estimates 
for the cold gas density and the SFR density are broadly 
consistent with the data. This supports our assumption that 
there exists a correlation between the declining cold gas mass 
with redshift and the declining efficiency of black hole for- 
mation with redshift. 



8 SIMPLE ESTIMATES OF THE QUASAR 
AND THE STAR FORMATION 
LUMINOSITY DENSITY IN THE UNIVERSE 

Extending the analysis of SFR density and cold gas density 
in the universe, we now compute the the relative contribu- 
tion to the ionising background arising from the star forma- 
tion rates and quasars. In the previous section we computed 
star formation rates by assuming a direct proportionality be- 
tween the fraction of cold gas in galaxies and the efficiency 
of black hole formation. In a more realistic scenario, the cold 
gas content in galaxies builds up gradually, and is controlled 
primarily by two competing processes. The first process is 
that of gas condensation, whereby the gas loses energy over 
a cooling time scale and falls into the central baryonic core 
over a dynamical infall time-scale to form stars. The sec- 
ond process involves the expulsion of a fraction of this gas 
from the galaxy as a result of feedback from supernovae and 
AGNs. It seems reasonable, therefore, to suppose that the 
gas content in galaxies evolves dynamically, peaking at some 
epoch and declining afterwards. 
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Figure 16. Left panel: The evolution of luminosity density (in 1500 A) for quasars with minimum black hole masses M,„i„ = 10^, 10** 
and 10^ Mq respectively. Right panel: The evolution of luminosity density (in 1500 A) estimated from the star formation rates in the 
universe (see text for details). The line-styles are same as in figure IT4l 



However, in our model the black hole formation effi- 
ciency does not peak at high redshifts but reaches a maxi- 
mum saturating value. Nevertheless, in view of the consis- 
tent results for the evolution of flgas and psfr{z), we bypass 
the exact evolution history of e{z) and fgas{z) by assuming 
their parallel evolution. This way we compute the maximal 
evolution of the luminosity densities for star formation rates 
and the quasars by extrapolating the model e{z) to high red- 
shifts. It is possible that at high redshifts, the evolution of 
star formation rates as well as the black hole formation rates 
are driven predominantly by a rap id gravitational growth of 
structure jPi Matteo et al1l2003l) . Therefore, at these red- 
shifts the exact evolution of e{z) and fgas{z) is perhaps not 
so important. 

Now the ultraviolet luminosity arising from star forma- 
tion can be written as t,Madau. Pozzetti fc Dickinson.,.1998il 

Luv = C X ( SFR/ Moyr"^) erg s-^Hz-\ (34) 

At 1500 A, just longward of the hydrogen Lyman limit (e.g. 
HNR98), the value of C for a Salpeter IMF is given as 
C ~ 8 X lO'^^. Thus the ultraviolet star formation luminosity 
density is written as 

/51500 = C X pSFR- (35) 

For quasars with a minimum black hole mass as Mmin, 
the luminosity density in B-band can be written as 

P4400{z) = dMBH X [ ^b{Mbh) ^3g^ 

X tgso{z) X Rbh{Mbh,z) ]. 

Here, tqso{z) x Rbh{Mbh,z) is the number of active 
quasars with black hole mass Mbh existing at redshift z. 



To convert the _B-band luminosity into luminosity at 1500A, 
we adopt a p ower law spectral energy distribution scaling 
as Li, oc v~°' jPeillQQ^l . where a — 0.5 and for a given 
band is simply Ltand/i^band- In figure ITCl we plot the evo- 
lution of luminosity densities due to star formation rates 
and quasars. These luminosity de nsities have also been est i- 
mated previously, for example, in lBovle fc TerlevichI l)l998() . 
For minimum black hole masses in the range lO'^ to 10^, we 
find that the drop in the quasar luminosi ty density towards 
low redshifts is steeper than computed in lBovle fc TerlevichI 
( 1998). For minimum quasar black hole masses of 10* Mq the 
curves in fig llfcil suggest that at high redshifts, the ionising 
background is predominantly of stellar origin. It comes out 
from our calculations that quasar dominated ionising pho- 
ton background is plausible if the minimum quasar black 
hole masses are < lO'^M© and the minimum circular ve- 
locities for haloes to effectively retain gas are fairly high 
(~ 150 km s~^). However, in light of the simplicity of the 
model used and the uncertainties involved in each of the 
estimates, these results should only be considered as indica- 
tive. 



9 DISCUSSION AND CONCLUSIONS 

We have presented an empirically motivated model for ac- 
cretion dominated self-regulated growth of super massive 
black holes in galaxies and investigated their implications 
for the evolution of the quasar population in the universe. 
The model looks into the core aspects of quasar population 
such as the space density evolution of quasars, the evolution 
of the quasar characteristic luminosity, plausible minimum 
masses of quasars, the mass function of the super massive 
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black holes and the formation epoch distribution of these 
black holes. 

The model suggests that the characteristic luminosity 
in the quasar luminosity functions arises primarily as a con- 
sequence of a characteristic mass-scale above which there is 
a systematic deviation of the black hole merging rates from 
the halo merger rates. Up to this mass-scale, the black hole 
merging appears to be closely following the merging of cor- 
responding dark matter haloes. This anomalous behaviour 
is probably a consequence of rarer galaxy merging in mas- 
sive haloes. Galaxies in massive haloes tend to retain their 
properties rather than quickly sink to the centre through 
dynamical friction. The extent of this effect might depend 
on the local environment of the merging galaxies. It is pos- 
sible, for example, that galactic merging in massive haloes 
is efficient for high sigma peaks only. We have mimicked 
such a behaviour by introducing a simple phenomenological 
decline in the black hole formation rate through a function 
Ffj,{M) with a power law drop at the halo high mass end. 
When combined with a declining efflciency of black hole for- 
mation with redshift which we attribute to a decrease in the 
cold gas mass available, the model can reproduce the quasar 
luminosity function for a wide range of redshifts. Further- 
more, we have shown that this declining black hole forma- 
tion efficiency gives a better comparison of the model break 
luminosity evolution with the evolution of characteristic lu- 
minosity used in the empirical double power-law fit to the 
data. 

On the lines of WL03, we have modelled time-scales 
for quasars proportional to the local dynamical time of the 
newly forming baryonic core, which in turn is proportional 
to the halo dynamical time. We found that the local dy- 
namical time gives a better comparison with low redshift 
observations than the Salpeter time-scale. Furthermore, the 
adopted redshift varying time-scale allows a natural trans- 
lation of the black hole formation rate density into number 
density of quasars, consistent with the the observed data. 

Calculation of the number density evolution of quasars 
in our model also led us to infer plausible minimum masses 
of optical quasars on the order of ~ 10* M©, consistent 
with the estimates of iSoltan il9fi3) "^ho computed minimum 
masses assuming a radiative efficiency of optical quasars 
^rad ~ 0.1. Additionally, our model suggests a less steep 
decline of quasar number density at high redshifts than is 
generally assumed. 

We also calculated the evolution of the comoving den- 
sity of black holes with masses above 10* Mq and 3 x 10* 
and 10^ Mq. The possible reasons for anomalous non- 
monotonic evolution of pbh inferred from the black hole 
mass function Ngjj have been discussed. We find that for 
minimum black hole masses of 10^ Mq a density of at least 
pBH ~ lO'^MoMpc"^ has assembled by redshift of 2. The 
absolute upper limit for minimum SMBH masses of 10*Mo 
is pBH ~ IO^MqMpc^'^. These estimates are consistent 
with the estimates of local black hole density (Chokshi & 
Turner 1992; Salucci et al. 1999; Yu & Tremaine 2002). Our 
model suggests that a significant fraction of SMBH mass 
density was already in place by z ~ 2 to 3. Interestingly, 
these estimates are also consistent with the recent results of 
iDi Matteo et al.l J2003h . Moreover, the model SMBH mass 
functions at 2 < 2.5 are c onsistent with the lo cal mass func- 
tion of SMBH derived in lSalucci et"ai] (Il999h . 



Finally, the predicted ionising background at high red- 
shift appears to be predominantly of stellar origin. We find 
that a quasar dominated ionising background would need a 
lower limit of minimum SMBH masses along with a higher 
value of minimum circular velocities of haloes to effectively 
retain gas in them. Obscured low mass AGNs could thus 
raise the prospects of a quasar dominated ionising back- 
ground. However, it is also possible that low mass AGNs 
have relatively low radiative efficiencies (e.g. Yu & Tremaine 
2002) and thus may not be able to contribute effectively to 
the ionising background. 
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APPENDIX A: A NOTE ON THE HALO 
MAJOR MERGER RATES 

The halo major merger rates derived in section (2.1) are 
based on the approach of Kitayama & Suto (1996) whereby 
the major mergers are characterised as merger events in 
which an object at least doubles its mass. As pointed out in 
section 2, in this definition of major merger events, there is 
always a chance that the main progenitor has more than half 
the mass of the final halo, even though it is not supposed 
to. We discuss here how one can systematically account for 
this possibility and compute corrected major merger rates 
where no progenitor has more than half the mass. 

Consider a merger event in which an object of mass M 
is formed from a progenitor of mass Mi (Afi < M/2). The 
mass accreted by this progenitor is M — Mi = AAf (AM > 
M/2). This mass AM can come either as a combination of 
formed smaller mass objects rrii such that ~ ^-^^ 

or as a single formed object of mass AAf. In order to dis- 
tribute the mass AM (variance denoted as 5'a) into smaller 
mass objects of mass m (variance denoted as Sm) we use 
the Press-Schechter mass function (at time t correspon ding 
to gi ~ Wi ~ aj2 ). In the excursion set picture of Bond ct alJ 
il99lfl . the random walks of trajectories 6{Sm) for the sub- 
clumps m now begin at Sm = 5'a instead of Sm = 0. This 
simply leads to an offset in the first passage distribution 



function, which can be written as 

f{Sm\S^)dSm = :^JS~S^: 



rexp 



2(5™ — Sa) 



dSn 



(Al) 

The probability that none of the masses nii is greater than 
M/2 can be inferred from the 'fraction' of mass of AM ex- 
isting in objects below mass Af/2. Thus 



Pm[m < (M/2) j AM] = erf 



V2{Sh - Sa) 



(A2) 



where 5";, = S{M/2). Note that when 5'a Su this prob- 
ability is one, that is all the masses are less than M/2, 
as expected. Hence the corrected major merger rate is given 
as 



Rmm{M, t) 



Nps{M,t) 



'dS,^P„ 
at 



(A3) 



= Nps{M,t) 
where I{M, i) is given as 



dui 



I{M,t) 



dSi 



2^ [Si - 5(M)]3/2 



I{M,t) 



erf 



;(t) 



y^2{Sh - Sa) 
(A4) 

Unlike the expression for merger rates given in KS96, we find 
that 7(M, t) here is a function of time. For any given mass 
M, therefore, KS96's expression becomes more accurate as 
redshift increases (or for large lo) as Pm — > 1. Unfortunately 
the full expression for I{M,t) is not easy to evaluate. The 
difference becomes more pronounced at low redshifts, even 
though for massive haloes (> W^^Mq at z = 0) the dis- 
crepancy is too weak to affect the results presented in this 
paper. Hence for the purpose of deriving quasar luminosity 
functions and other properties we have restricted ourselves 
to the 'computationally easy' expression of KS96. 

An additional advantage of equation lAll is that it could 
be used in the merger trees to sub-divide the remaining mass 
after picking the main progenitor at an earlier time in a short 
time-step. If this time step is very small, then in most cases 
the remainder mass is small and can be treated as accreted. 
In other cases, we think it more satisfactory to distribute the 
remaining mass (other than the main progenitor) according 
to equation lAll rather than use the various other methods 
that have been proposed in the literature. 
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